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1. Introduction 

Pion-nucleon scattering is a fundamental reaction that gives access to fundamental questions 
related to strong interactions. At higher energies, for example, allow us to study the baryonic 
spectrum of QCD and its properties. On the other hand, at low energies it is an excellent test for 
the chiral dynamics of QCD , able to provide a systematic framework to study isospin violation as 
well as valuable information about the internal structure of the nucleon. Regarding the latter, there 
is a strong demand from the dark matter community for an accurate value of the scalar coupling 
of the nucleon at zero momentum transfer, which is the definition of the pion-nucleon sigma term 
(cTjtiv)- This is so because there are important discrepancies between the values for o k n obtained by 
different partial wave analyses (PWAs). Namely, the classical partial wave analysis (PWA) of the 
Karlsruhe group gives a value of g k n « 45(8) MeV while the more modern PWA of the George 
Washington University group obtains g k n ~ 64(7) MeV [Q]. This disagreement on the sigma-term 
value is the main hadronic uncertainty for the direct detection of dark matter [Qj. 

Chiral perturbation theory (ChPT) is an excellent tool that can shed light on this issue because 
it provides a systematic and model independent way to study perturbatively the strong interactions 
of hadrons. From this perturbative construction one can also include formally, in a quantum field 
approach, the interaction with scalar sources [Q, from which the scalar form factor of the nucleon 
can be investigated. In fact, ChPT has been used several times to estimate the value of the pion- 
nucleon sigma term (see Refs. ^ ^]), although the poor convergence of the chiral series prevented 
to make any conclusion about its value. As we show here, the lack of convergence of the previous 
analyses can be overcome if one works within a framework where one preserves the good analytical 
properties of a covariant calculation and includes the relevant degrees of freedom for the process 
considered. Namely, for nN scattering, it is known that the A(1232) resonance plays a fundamental 
role because its closeness to the nN threshold. This makes that the behaviour of this resonance 
cannot be faithfully reproduced by a finite polynomial via the resonance saturation hypothesis [ph. 
In our work we considered explicitly the contributions of the A(1232) following [§]. 

One added difficulty that one should overcome when working with covariant baryon chiral 
perturbation theory is the breaking of the perturbative expansion already addressed in Ref. 
This problem is solved when working within the extended-on-mass-shell scheme (EOMS), which 
cancels the PCBT (which are analytical pieces) by renormalizing the low-energy counterterms 
(LECs) of the most general chiral Lagrangian [ |To| , [FT]]. 

One of the main advantages of this scheme respect to the other covariant approach, infrared 
regularization (IR) JT2[], is that EOMS preserves the good analytical properties of the relativistic 
scattering amplitudes. However, IR gives rise to unphysical cuts that limit the convergence of the 



chiral amplitude, as shown in [ ]12| , |13| , |14| ]. For this reason, we decided to work within the EOMS 
to preserve the good analytical properties of our calculated amplitudes as well as the standard 
power counting of ChPT. The results presented here are part of the detailed analysis of the nN 



scattering phenomenology of Ref. Q14|], where the nN scattering amplitude is calculated up to 
- 3 ) including the A( 1232). 
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2. Fits 



In order to determine the LECs, which numerical values are not fixed by chiral symmetry, we 
decided to fit our chiral amplitude to the phase shifts provided by three different PWAs. These are 
the PWAs of the Karlsruhe group (KA85) [jl|], the George Washington University group (WI08) [§] 



and the Matsinos' group (EM06) [[16J]. 

In Fig. p] we observe that in the 4-ChPT case (green dashed line) we are able to fit well the 
WI08 phase shifts up to energies of y/s = 1.14-1.16 GeV, depending of which partial wave one 
consider. This situation is considerably improved once the A(1232) is included explicitly in our 
formalism A-ChPT (red solid line), in which case we are able to describe the data very well up to 
energies of y/s = 1.20 GeV. Although not shown here, this good agreement is also observed for 



KA85 and EM06 (see Ref. [14]) 
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Figure 1: Fits to WI08 [g] with (solid red line) and without (green dashed line) the inclusion of the A(1232) 
as an explicit degree of freedom. For 



The values of the LECs extracted from these fits are shown in Table [j] The results of the 
A^ChPT are given in columns 5-7, while the results of the A-ChPT fits are shown in columns 2-4. 
In the latter case we have one extra LEC compared to the ^-ChPT, the NA axial coupling Qia). 
This LEC, on the other hand, can be related to the A(1232) Breit-Wigner width, So, taking as 
input the value reported by the PDG T A = 1 18(2) MeV 0], one deduces a value of h A = 2.90(2). 
However, we decided to leave this LEC as a free parameter to check the reliability of the different 
PWAs. Once we have determined completely the chiral amplitude, a key point is to predict other 
interesting phenomenology. Here we present the results regarding the convergence of the chiral 
amplitude in the subthreshold region, Sec. || and the extraction of g k n, Sec. || 
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0.11 
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Table 1: This table gathers the values of the LECs extracted in the different fits, with and without the 
A(1232). The @(p 2 ) and ^(/? 3 ) LECs (the a and d t ) are shown in units of GeV 1 and GeV 2 respectively. 



3. Subthreshold region 



The subthreshold region is of special interest because contains points that are related to impor- 
tant low energy theorems. One of the most important and known of these theorems in nN scattering 



is the Cheng-Dashen theorem [18], that relates the value of the Born-subtracted isoscalar scatter- 
ing amplitude, D + , at the Cheng-Dashen point, to the scalar form factor of the nucleon, o{t), at 
t = 2M\. This theorem is used by several PWAs to extract the value of o n s, since the difference 
a(2M 2 ) — Ojzn = A a can be derived in different ways, see for example Refs. [^, |5[ 19]. A dif- 
ficulty one had in the subthreshold region is that the ChPT analyses of nN scattering could not 
extract, from physical information, the value of the subthreshold quantities that the PWAs are able 
to obtain by means of dispersive methods. Therefore a good way to explore the convergence of 
the chiral series in this region is to compare our results there with the one obtained by their corre- 
sponding PWAs. The subthreshold quantities that we studied are the subthreshold coefficients d^ Q 

Jqq + d^t + . . . and the so called E-term, 



and cIq V defined by the power expansion D + (v = 0,f) 

£ = f^D + (v = 0,t = 2M 2 .). We chose these quantities because they are closely connected to the 



value of Ojin [14] 
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1.73(19) 
84(10)' 


-1.65(28) 
1.70(18) 
103(5)' 


-1.56(5) 
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Table 2: Values of the subthreshold quantities under consideration, extracted from direct extrapolation of 
the chiral amplitude into the subthreshold region. They are compared with the values reported by the PWAs 
(last two columns). 



In Table ^ we show the result of the extractions of d^ , d^ and £. We show there that, in the A- 
less case, the general trend is to overestimate, in modulus, the value of these subthreshold quantities 
compared with those calculated by the corresponding PWAs. However, including the A as an 
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explicit degree of freedom, we are able to extract values of the subthreshold coefficients and the E- 
term that are in agreement with the PWAs. This suggest that the A(1232) is a key ingredient for the 
convergence of the chiral series also in the subthreshold region. Once this resonance is included as 
an explicit degree of freedom we could, for the first time in the literature, connect the information 
that lies in the physical region, encoded in the LECs through our fits, with the amplitudes in the 
subthreshold region. This information was, up to now, only accessible by dispersive methods. 
However, the A(1232) alone is not enough to achieve the best convergence, and is also necessary 
to work in a formalism where we preserve the good analytical properties of a covariant calculation 
to extract all the potential of baryon ChPT (see Ref. JH] ] for more details). 



4. The pion-nucleon sigma term 

The pion-nucleon sigma term is an important quantity that contains information about the 
internal scalar structure of the nucleon. It is related to the origin of the mass of ordinary matter and 
is important for the investigation of QCD phase diagram. In addition, nowadays it is widely used 
for estimations of the dark matter-nucleon spin independent elastic scattering cross section, which 
are used for direct detection of dark matter. 

The PWAs can extract the value of g k n by extrapolating the D + to the Cheng-Dashen point 
(v = 0,? = 2M^). 2 However, this is a delicate extrapolation into the subthreshold region with no 
direct physical information to compare with. In this regard, the advantage of ChPT over dispersive 
methods is that chiral symmetry allow us to relate g k n to the LEC c\ (in an 0(p 3 ) calculation) 
which, on the other hand, can be determined from experimental information. That means that once 
the value of ci is reliably fixed employing experimental information, and the convergence of the 
chiral expansion for g„n is proven, one may give a reliable extraction of G n ^- 

The explicit relation between g„n and c\ can be obtained directly form a(t = 0) or by applying 
the Hellmann-Feynman theorem on the chiral expansion of the nucleon mass. By both means one 
obtains, in an & (p 3 ) covariant calculation, the result 3 [[Z(]] : 




G KN = -A Cl Ml- //f^ \ arccos^+M.log^ | (4.1) 

The point here is that our good convergence of the chiral series above and below threshold, 
once the A(1232) is included, give us confidence about the reliability of the LECs extracted with 
A-ChPT from PWA phase shifts. On the other hand, the chiral expansion for g k n exhibits a good 
convergence up to the order that we calculate, as was already shown in Ref. [p0[|. Using the values 
of C[ extracted with A-ChPT, one obtains the values for g^n shown in Table |3[ 



'These values of £ are extracted from an ff(p^) ChPT calculation, which is known that underestimates the value 
of this quantity in approximately 10 MeV |l9(]. As discussed in Ref.[14] what matters for determining c%;v is = 
fi(d++2Mld+). 

2 The variable v is defined in terms of the Mandelstam variables s and u as v = 



3 In Eq. ( t.l ), cj corresponds to the EOMS renormalized LEC. 
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KA85 
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[§ 


EM06 
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ci (GeV- 1 ) 


-0.80(6) 


-1.00(4) 


-1.00(1) 








o KN (MeV) 


43(5) 


59(4) 


59(2) 


45(8) 


64(7) 


56(9) 



Table 3: Values of g k n extracted from the fits to the different PWAs. 



As one can see in Table |3[ the A-ChPT extraction is in good agreement with the one of the 
PWAs. This suggests that the discrepancy between KA85 and WI08 regarding the value of g^n is 
not due the methodology used by these PWAs. Moreover, we observe that our extractions employ- 
ing the modern PWAs WI08 and EM06 agree remarkably well, despite both follow a veiy different 
systematics. However, what WI08 and EM06 have clearly in common is that both employ modern 
and high quality data. This suggests that the modern data point to a relatively large value of G % ^. 

In order to discriminate which extractions are more reliable we decided to calculate different 
quantities that can be compared with independent determinations. Two of those quantities are h& 
and the Goldberger-Treiman deviation, Act- For Ha we showed that only using the WI08 solution 
to fix the LECs, we extract a value for this coupling that is related to a Ta compatible with the 
value reported by the PDG. In the case of A^r, our extractions are compatible with independent 
determinations from AW scattering and pionc-atoms only if we employ the modern PWAs of WI08 
and EM06. Nevertheless, the most important quantity to compare with is maybe the scalar-isoscalar 
scattering length (%,), since this quantity is closely related to the value of g k n- In order to illustrate 
this, we show in Fig. ||] the dependence of £^ (Lj — g k n ~ 3 MeV) as a function of and a\ + . 
This relation was already shown in Ref. J2lj|, from where Gasser, Leutwyler and Sainio deduced 
a value for g k n of 45 MeV. In fact, the value of Oq + was key in their extraction (see Fig. 1 of 
that reference). There is clearly shown that the value of g k n that they report (point A) relies on a 
negative value of a^ + , accepted by the time in which this paper was published. However, according 
to modem determinations for this scattering length extracted from pionic-atom [23], the value of 
a~Q + is now well consistent with zero. This affects directly to the value of g k n, which is pushed to 
larger values. This means that using exactly the same argument that Gasser, Leutwyler and Sainio 
employed to determine a value of g k n ~ 45 MeV would lead to a larger value of the sigma term in 
light of this updated determination. 

We also used the value of a^ + extracted from 7r-atoms data to check the reliability of the PWAs 
employed in our extractions. As we show in Table fj], we observe, again, that only using modern 
PWAs our extractions are more compatible with independent determinations. 

From these checks with independent and updated phenomenology, we conclude that the most 
reliable PWAs to use as input to extract g^n are the modern PWAs of the George Washington Uni- 
versity and Matsinos' groups. From these modern PWAs we extract a value of g k n = 59(7) MeV, 
where the 7 MeV of error takes into account the uncertainty in the LEC c\ and the first higher order 
correction to g k n not taken into account in our 0(p 3 ) calculation. 

This value, unlike to the old value of 45 MeV, is consistent with the updated phenomenology 
related to the sigma term (see Ref, [[14], 20, 24]). This consistency with updated phenomenology is 
important to claim reliably for a certain value of g k n and, unfortunately, is not satisfied by other 
experimental based determinations [25]. The more recent calculation of Ref. [26], lacks also from 
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consistency with independent determinations 5 and extract results at odds with the PWA that they 
employ. 




+0.01 0.00 -0.01 -0.02 



4+ [AC 1 ] 

Figure 2: Dependence of E f / (diagonal lines) with with a^, and a\ + . The point A correspond to the values 
of Oq + and a\ + deduced from the data of Bertin et. al. p2| ] with one standard deviation (ellipse). The point 
C and the square correspond to the central values and errors for a^ + and af, reported by KA85. 





KA85 A-ChPT[14 


I 


WI08 A-ChPT[14 


] 


EM06 A-ChPT ^4 


] 


7T-atoms jj23] 


a+ (10- 3 M^) 


-11(10) 


-1.2(3.3) 


2.3(2.0) 


-1.0(9) 



Table 4: Values of aj, extracted from the different PWAs (columns 2-4) compared with the independent 
determination of Ref. p3|]. The value shown in column 4 is obtained from Ref. [^3j taking only n + p and 
K~ p scattering data. 



5. Summary and Conclusions 

In this contribution we highlighted some of the results obtained in Ref. [|l4|]. There, we cal- 
culated the TtN scattering amplitude up to &(p 3 ) in the EOMS covariant scheme including the 
A(1232)-resonance as a dynamical degree of freedom. We show how with this inclusion one 
achieves the best convergence in the physical and the subthreshold regions. Thanks to this im- 
provement we were able to connect, for first time in the literature, the information extracted in the 
physical region with the one that lies in the subthreshold region. This good convergence of the 
the chiral amplitude in both regions allowed us to extract reliably the value of the pion-nucleon 
sigma term form experimental information (PWAs), since the value of o k n converges well in our 
covariant &(p 3 ) calculation including the A(1232) [[14], ^ E7|]. Employing the modern PWAs we 
extracted a value of g„n = 59(7) MeV, and compared this value with related phenomenology. We 
show how this relatively large value of o k n is consistent with the updated phenomenology, while 
the old value of 45 MeV is only constent with a negative value of a^ + , which is at odds with mod- 
ern 7T-atoms results. This consistency with updated phenomenology gives a strong support to the 

5 Namely, the Ha extracted in this analysis is not compatible with the result deduced from the PDG. The threshold 
parameters, on the other hand, are not even considered in this work. 
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relatively large values of O n x, that has been recently shown to be also consistent with a negligible 



strange quark contribution to the nucleon Q24Q. These updated determinations of the sigma terms 



(Ref. [|(| Hi) from effective field theory can be very useful to the dark matter community. 
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